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We describe a method for improving coherent control through the 
use of detailed knowledge of the system's Hamiltonian. Precise uni- 
tary transformations were obtained by strongly modulating the sys- 
tem's dynamics to average out unwanted evolution. With the aid of 
numerical search methods, pulsed irradiation schemes are obtained 
that perform accurate, arbitrary, selective gates on multi-qubit sys- 
tems. Compared to low power selective pulses, which cannot aver- 
age out all unwanted evolution, these pulses are substantially shorter 
in time, thereby reducing the effects of relaxation. Liquid- state NMR 
techniques on homonuclear spin systems are used to demonstrate the 
accuracy of these gates both in simulation and experiment. Simu- 
lations of the coherent evolution of a 3-qubit system show that the 
control sequences faithfully implement the unitary operations, typi- 
cally yielding gate fidelities on the order of 0.999 and, for some se- 
quences, up to 0.9997. The experimentally determined density ma- 
trices resulting from the application of different control sequences on 
a 3-spin system have overlaps of up to 0.99 with the expected states, 
confirming the quality of the experimental implementation. 



I. INTRODUCTION 

The past decade has seen a substantial interest in improv- 
ing coherent quantum control. Coherent control had origins 
in both nuclear magnetic resonance (NMR) [jl],^ and opti- 
cal spectroscopy [|3jf5J] . For an overview of advances in both 
fields, the reader is referred to [|5|] . Since coherent control's in- 
ception, many different techniques have been used both to im- 
prove selectivity and to reduce the duration of control pulses. 
For spin systems, the Fourier transform has been used to ap- 
proximate the excitation profile in the limit of low power and 
no spin-spin couplings |g] and more complete analytic so- 
lutions have been developed to aid in general pulse design 
and analysis [f7|-|l2|] . Alternatively, very sophisticated shaped 
pulses have been designed using a variety of computer-aided 
methods [ |l3|]l4| ] or feedback from system observables Jl5|]. 
Equivalent analytic theories jT^ , [l7|], computer-aided meth- 
ods [ |l8| ] and feedback techniques [ |19p0| ] have also been de- 
veloped by the optics community. Similar techniques are also 
used in other fields such as the control of trapped ions 

The development of liquid- state NMR systems as prototype 
quantum information processors [ ^3| , ^4| ] has enabled experi- 
mental demonstration of quantum algorithms [[Z5|-|27|], quan- 
tum error correction [ |28[|2^ ], and quantum simulations j3Q|] . 
These experiments built upon well-established spectroscopic 
techniques developed over the past four decades, such as using 
low-power (soft) shaped radio-frequency (RF) pulses to obtain 
selective operations. However, the selective pulses employed 



to date have the disadvantage that low power implies long du- 
ration. This not only introduces errors due to relaxation, or 
decoherence, but also allows significant evolution under the 
action of the internal Hamiltonian. In the past, this evolu- 
tion was rarely of concern because there was little importance 
placed on implementing a particular operation. For example, 
in spectroscopy there are entire classes of propagators that se- 
lectively excite a single spin from its equilibrium state, but for 
applications such as quantum computing the transformation 
must act as expected for all input states. A second problem 
with soft-pulse techniques is that selective pulses simultane- 
ously applied to different spins interfere with each other, thus 
causing significant deviation from the desired action [31]. To 
address some of these problems, several groups precalculate 
these errors and incorporate corrections into their analysis and 
pulse design [ f32[|33| ] . However, not all errors can be corrected 
using these techniques [[34|], and it would be preferable to aver- 
age out unwanted evolution by the use of strong control fields, 
so that no additional corrections are required. 

In this paper, we present a procedure for finding high- 
power pulses that strongly modulate the system's dynamics 
to produce precisely a desired spin- selective unitary propa- 
gator. These operations, or gates, allow arbitrary rotations of 
each spin around independent single- spin axes, while refocus- 
ing the internal evolution. They are "self-contained," in the 
sense that they can be placed back to back in longer sequences 
without requiring additional computational resources or post- 
experiment corrections. By using high-power, pulse durations 
are decreased by almost an order of magnitude, thereby sig- 
nificantly reducing the effects of relaxation. Finally, the use 
of strong modulation also opens the possibility of incorpo- 
rating robustness against slowly varying or time independent 
incoherent errors such as those caused by static field or RF 
inhomogeneities [[35|-[37|] . Our control methods are the first to 
combine all of these features. 

The pulses presented here have been applied in re- 
cent Quantum Information Processing (QIP) experiments 
to demonstrate algorithms [[27]], study notions of measure- 
ment [|38|], and test new methods for noise control [39]. They 
promise to be increasingly useful in future NMR QIP exper- 
iments, where larger numbers of qubits will necessitate in- 
creasing the number of homonuclear spins. In addition, these 
methods can be adapted to develop improved pulses for se- 
lective spectroscopy and imaging pl| ] . Finally, although 
presented within the context of NMR, these methods are ap- 
plicable to any system where the total Hamiltonian is well 
known and the external degrees of freedom allow for universal 
control, both requirements of any quantum information pro- 
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cessor. 



II. CONTROL OF QUANTUM INFORMATION 

In the standard model of quantum computing, an algorithm 
can be expressed as a series of unitary operations that maps 
a set of input states to a particular set of output states. The 
physical implementation of an algorithm requires the use of a 
quantum system with a Hamiltonian that contains a sufficient 
set of externally controlled parameters to allow for the gener- 
ation of a universal set of gates [42]. The task of control is 
to find a time-dependent sequence of values for these control 
parameters that modulates the system's dynamics in order to 
generate a particular gate to the required precision. 

Given a control sequence, solving for the effective Hamil- 
tonian is straightforward. Unfortunately, going the other way 
is much more difficult. That is, finding an RF waveform that 
produces a propagator with desired properties is an inverse 
problem. Traditionally, analytic techniques, such as average 
Hamiltonian theory [g], have been used to determine an ap- 
propriate control sequence. With modern computer resources, 
numerical methods provide a more efficient and accurate so- 
lution to this problem. 



A. Gate Fidelity as a Metric for Control 

A metric of a gate's performance should describe the qual- 
ity of a general transformation, including the possibility of 
non-unitary evolution. Unfortunately, such information is not 
conveniently accessible by experiment, so we choose a metric 
comprised only of sets of state measurements. For an input 
state, pi n , the ideal transformation maps the system to a theo- 
retical output state, p t h, i-e., 



Pir, 



Pth. 



On the other hand, a simulated or experimentally implemented 
control sequence will produce a different output state, p out , 
i.e., 



Pin 



Pout- 



(2) 



Noting that p is Hermitian, the projection between these two 
states, defined as 



P(Pth,Pout) 



trace(p th p out ) 



y/trace(p 2 th )trace(p 2 out ) ' 



(3) 



quantifies how similar in 'direction' the two states are. This 
metric is analogous to the dot product between two vectors, 
varying from —1 for anti-parallel states to 1 for identical 
states. A value of zero indicates orthogonal density matrices. 
In order to account for non-unitary evolution, a second term 
multiplies the projection yielding the attenuated correlation, 
namely, 



C(pth, Pout) = P(pth, Pout) 



I trace(p 2 out ) 



trace(p 2 n ) 
_ trace(p th p out ) 
yjtrace(p 2 th )trace(p 2 n ) 

We define the gate fidelity, F, of a transformation as 



F = C(pth,Pout), 



(4) 
(5) 

(6) 



where C represents the average attenuated correlation 
over an orthonormal set of input density operators (i.e., 
Trace[pjpk) = Sjk) that span the Hilbert space. It should 
be noted that F is maximized (with a value of one) when the 
implemented and ideal transformations are the same, and is in- 
sensitive to differences in the global phase between the ideal 
and implemented transformation. 

We can derive a useful alternate form for the gate fidelity 
in terms of the actual and theoretical transformations instead 
of the input-output state relations. This form is both eas- 
ier to compute and has intuitive appeal in that knowledge of 
the transformation can be directly translated to gate fideli- 
ties. First we assume that our ideal transformation is unitary, 
and the implemented transformation is a completely positive, 
trace-preserving linear map . In other words, the imple- 
mented transformation takes normalized density operators to 
normalized density operators and if the system starts as sub- 
system of an entangled system, then the full system's density 
operator also maps in a reasonable way. Under these assump- 
tions, Eq. ([]]) and (||) are explicitly given by 



Pth = UthPinU. 



thi 



and 



Pout 



(!) where the satisfy 



£4^ = i- 



(7) 



(8) 



(9) 



We now show that the gate fidelity reduces to 

F = J2 \Trace(UlA^)/N\ 2 , (10) 

/i 

where N is the dimension of the Hilbert space. Using the nor- 
malized Pauli basis, aj, as the orthonormal input density op- 
erators and the cyclic properties of the trace, Eq. (|^) becomes 

N 2 

F = J2Tr[(U th a 3 Ul)(^2 A^A^/N 2 (11) 

N 2 

= £ Tr[aj £ UlA^jAlUtnyN 2 . (12) 

Expanding the product of U^A^ in terms of the orthonormal 
Pauli basis (U} h A„ = £ fe B*a k ), yields 
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= Y, Tr fo(E B >^jCE B ™*°rn)]/N 2 (13) 

jfi k m 

= B k ^*Tr[a 3 a k a 3 a m }/N 2 . (14) 

j fikm 

Because the a basis is orthogonal, only terms where k = m 
contribute. Therefore, Eq. ( |l4| ) reduces to 



F = Y J \Bl\ 2 Tr[cj J a k a 3 a k \/N 2 

jfik 



(15) 



If is not proportional to identity, it will anti-commute with 
exactly half the <jj terms in the sum, while commuting with 
the other half. Therefore, two sets of terms cancel and have 
no contribution to F. Defining g\ to be the element that is pro- 
portional to identity, Eq. ( ^5[ ) further reduces to 



F = Y,K\ 2 Tr[o 3 a 3 \/N* 



Y^\bI\ 2 /n^ = y.K\ 2 / n - 

3fJ> M 



(16) 
(17) 



This is clearly equal to Eq. [U]. Thus, the gate fidelity cor- 
responds to how well the actual transformation reverses the 
action of u\ h . In this form it is obvious that the gate Fidelity 
is independent of which orthonormal basis of input states are 
used as pi n . 



B. NMR as an Example System 

As an example, liquid-state NMR is used to demonstrate 
how to find control sequences to implement particular gates. 
In NMR, spins in a large static magnetic field (in our case, 
9.7 T) are controlled via external RF pulses. The internal 
spin Hamiltonian is composed of both Zeeman interactions 
with the applied field modified by electron screening (chem- 
ical shift) and scalar couplings with other spins. Together 
these provide the QIP requirements of addressability and con- 
ditional logic respectively. In terms of spin operators, the in- 
ternal Hamiltonian is 



H in t = + 2^^ J kj I k 'I j , 



(18) 



k=i 



j>k k=l 



where uok represent the chemical shifts of the spins, Jkj the 
coupling constant between spins k and j, and n is the number 
of spins. The external Hamiltonian describing the coupling 
between the spins and an oscillating RF field generated by a 
single transmitter is 

n 



k = l 



(19) 



where ujrf is the transmitter's angular frequency, (j) the initial 
phase, and uj the power [p4|]. Of course, additional species 



can be added by including appropriate terms in H int and an 
additional H ext for each additional RF field. 

Using this knowledge of the internal Hamiltonian and the 
form of the external Hamiltonian, the parameter values that 
generate the desired gate must be determined. Here, a quality 
factor Q = 1 — \ /r F is minimized by searching through the 
mathematical parameter space using the Nelder-Mead Sim- 
plex algorithm [45]. While this function has many local min- 
ima, the Simplex algorithm often succeeds in finding satis- 
factory solutions. Our goal is to show that sufficient, imple- 
mentable control sequences can be found. Finding the opti- 
mal solution is much more challenging and based on our sys- 
tem and control parameter values, is not expected to improve 
pulse performance significantly. We have parameterized the 
control sequence as a cascade of RF pulses with fixed power, 
transmitter frequency, initial phase, and pulse duration (r). As 
will be seen, this is a particularly convenient and completely 
general parameterization, but we make no claims that it is the 
only, nor necessarily the best choice. If the RF power is con- 
stant over the duration of a pulse, i.e., the pulse's amplitude 
is square, the total Hamiltonian H tot = H int + H ext can be 
made time independent by transforming into the frame that 
rotates at the frequency of the transmitter. This allows the 
Liouville-von Neumann equation of motion to be solved by a 
single diagonalization. Initially, the starting density matrix is 
the same in both frames (p(0) = p(0)), so that at the end of 
the pulse, the density matrix in the new frame is given by 



p(r) = e 



-iH e 



fT p(0)e 



(20) 



where H e ff is the effective Hamiltonian in the new 
frame [ fl6| ] . Transforming this density matrix back to the orig- 
inal rotating frame gives 



p = U z (T)- 1 e- iH °" T p(0)e tIi °" T U z (T), 



iH eff ri 



where 



U z (r) = (e* 



)• 



(21) 



(22) 



Therefore, in the original rotating frame, the transformation is 
given by 



U period (r) =U- L (r)e- iH *ff r . 



(23) 



Because the evolution under the whole sequence is given in 
the original rotating frame, no additional resources are re- 
quired to concatenate pulses, nor is any mathematical correc- 
tion required at the end of an experiment. 

Cascading these periods yields the net transformation 



N 



U net = H U-\r m )e- iH ^ m ^^ m ^ 



(24) 



where the index m refers to the m th period, i.e., to the m th 
square pulse, with a corresponding set of 4 parameters. In 
other words, N constant amplitude pulse periods, each with 
a different transmitter frequency and initial phase, are applied 
in series. Clearly, a single period is not sufficient to generate 
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an arbitrary transformation; therefore the number of periods is 
increased until a suitable net transformation is found. Using 
desktop computing resources, this yields convergence times 
for three- and four-spin systems that are typically seconds to 
minutes. 

In addition to the desired propagator, JJ r ideal, an initial set 
of starting parameters for the pulse shape is required. While 
this initial guess must be reasonable {i.e., in the vicinity of the 
solution), many different starting points typically converge to 
equally deep minima. We have observed that the number of 
acceptable solutions for this parameterization is very large, al- 
lowing experimental implementation issues to be considered. 
For example, experimental limitations do not allow arbitrarily 
high powers or frequencies to be implemented. To keep the al- 
gorithm from returning infeasible solutions, a penalty function 
that increases as the parameter value moves towards infeasible 
solutions is added to the quality factor. Penalty functions are 
also used to guide the algorithm towards more favorable pulse 
solutions. In our case, penalties are placed on high powers, 
large frequencies, and negative- or long-time periods. 



III. SIMULATIONS 

The methods described above were used to obtain a set of 
pulses that implement each of a set of important single- spin 
gates. To study the performance of these gates, propagators 
for each of the pulses were simulated under different condi- 
tions. First, the performance of the pulses under idealized 
experimental conditions is considered. Second, the gate fi- 
delity is simulated as a function of systematic distortions of 
the pulse parameters. From these results, the relative impor- 
tance of implementation precision is determined. Finally, sim- 
ulations show that a pulse generates quite different evolutions 
as a test spin's resonant frequency is varied over a range of 
chemical shifts. 



A. Ideal Pulse Simulations 



needed vary from 0.9999 to 0.999999 depending on the as- 
sumptions used. These simulations predict an achievable level 
of control that approaches the most optimistic estimates for 
fault tolerant computation. As expected, the pulse duration 
decreases with increasing chemical shifts dispersion (selec- 
tivity condition) and, for the case that Jjk « \^k — \, the 
fidelity of the sequence decreases with increasing ratio of the 
couplings (bilinear terms) to chemical shift. 



B. Variations in the External Hamiltonian 

The external RF parameters are determined by a minimiza- 
tion procedure, suggesting that small variations of the exter- 
nal parameters should have little effect on the quality of the 
pulses. To check this assumption, the gate fidelity was calcu- 
lated as each of the six pairs of the four control parameters 
were varied over a range of errors typical of an experimental 
implementation. As a sample set, one pulse for each of the 
two systems is presented here. The results shown in Fig. || 
demonstrate the natural robustness against typical variations 
in the initial phase, frequency, and duration of each period. 
Clearly, the sequence is most sensitive to power variations. 
For the pulses listed in Table [| if the power's amplitude is 
changed by 5% the average fidelity falls to 0.96 ± 0.01 for 
Alanine pulses and 0.94 ± 0.04 for Crotonic acid pulses. For 
the 25 pulses used in [ |3^ ] the average fidelity at 5% amplitude 
deviation is 0.97 ± 0.02. For 10% deviation, the gate fidelity 
drops to 0.86 ±0.03 for the Alanine pulses and 0.81 ±0.12 for 
the Crotonic acid pulses. This pulse sensitivity to RF ampli- 
tude suggests that RF inhomogeneity may be a leading cause 
of experimental errors. While techniques to select homoge- 
neous regions are available [ [33|j48| ], the loss in signal to noise 
is significant, especially if multiple coils are used. Instead, 
because these errors are incoherent in nature, it is possible to 
design pulse sequences that refocus such homogeneities. Pre- 
liminary results confirm that sequences can be re-optimized in 
the presence of significant inhomogeneity. 



Pulses were created for three- ( 13 C-labeled Alanine) and 
four- ( 13 C-labeled Crotonic acid) spin homonuclear systems. 
The chemical shifts and scalar coupling constants for each of 
these systems are listed in Fig. [j]. As a representative set, each 
of the single spin tt/2, and nearest-neighbor paired tt pulses 
were simulated with the relevant characteristics summarized 
in Table | and example waveforms shown in Fig. ^[ The dura- 
tion of the pulses are on the order of 200/iS for the three-spin 
system and 420/iS for the four- spin system, both significantly 
shorter than those that could be obtained using low-power 
pulses. The average fidelities for each system are 0.9995 and 
0.995, demonstrating that, at least under ideal conditions, con- 
trol sequences that implement the desired transformation with 
high fidelity can be found. 

The ultimate goal of control in quantum computing is to at- 
tain fault- tolerant computation. While it has been proven that 
perfect control is not required [|47|], estimates of the precision 



C. Variations in the Internal Hamiltonian 

For NMR spectroscopy, the goal is to excite selectively 
spins in a band of frequencies leaving all other possible spins 
(with unknown precession frequencies) along the z axis. This 
requires that the propagator for spins at any other frequency 
be, at most a phase change. With detailed knowledge of the 
internal Hamiltonian, the effect of the applied RF field needs 
only be considered at the resonance frequencies of the chemi- 
cal species present in the given molecule. By relaxing the re- 
quirement that the effective Hamiltonian be zero for all chem- 
ical shifts other than those in the band of excitation, an RF 
shape can be found that more efficiently implements the de- 
sired gate for the frequencies of concern yielding high-power 
yet selective pulses. To demonstrate this idea more clearly, 
the gate fidelities of the two sample pulses considered in the 
sub- section B were calculated as a function of a test spin's 
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resonant frequency. As can be seen in Fig. [|, the fidelity is 
close to unity only near the resonance frequency for which 
the pulse was designed to work. This stresses the necessity 
of having accurate knowledge of the system's Hamiltonian. 
On the other hand, looking at the region immediately around 
the resonance we see the fidelity falls off quite slowly. This 
implies that small variations in the chemical shift do not sig- 
nificantly affect the fidelity of the pulses. For example, in 
the experiments presented below, the unwanted scalar cou- 
plings to the hydrogen atoms, which are equivalent to errors 
in the resonance frequency, were automatically refocused by 
the control pulse. It should be noted that no constraint was 
used to require this robustness, but that it results from the use 
of strongly modulating pulses. If this natural robustness is not 
sufficient, additional constraints can be added. Of course, this 
robustness also implies that selective pulses will be harder to 
design for spin systems where the spectrum is dense. 



IV. EXPERIMENTAL DEMONSTRATIONS 

In this section we experimentally demonstrate the efficacy 
of a representative set of pulses via both spectra and recon- 
structed density matrices. Each waveform was discretized by 
sampling at a constant rate that was faster than the largest 
phase modulation during any of the periods (sample times 
of order 0.5/is). Simulations of the sampled waveform con- 
firmed that the decrease in the gate fidelity was acceptably 
small. All experimental tests were carried out on the three 
carbons of 13 C-labeled Alanine (see Fig. |]for internal Hamil- 
tonian parameters) using a Bruker Avance spectrometer. 

First, a series of selective single- spin pulses were applied to 
the thermal state to demonstrate that these selective rotations 
fully refocus the internal Hamiltonian and so can be concate- 
nated arbitrarily. Fig. || shows sample sequences and the re- 
sulting spectra. In addition, because the internal Hamiltonian 
is fully refocused, applying selective transformations on dif- 
ferent spins sequentially has an effect equivalent to applying 
all the appropriate transformations to each of the spins simul- 
taneously (neglecting relaxation). With simultaneous, fully 
self-refocused tt pulses, selective couplings can be efficiently 
implemented using previously published techniques [^6 49], 
and is demonstrated in Fig. ^[ 

Second, the projection and attenuated correlation between 
the expected and experimentally determined results of differ- 
ent control sequences were measured for a set of input states. 
While spectra contain information about the observables of 
the current state of the system (single-spin transitions in the 
case of NMR), a single spectrum does not contain enough in- 
formation to reconstruct the entire state of the system (density 
matrix). By using "read-out" pulses that rotate different el- 
ements of the density matrix into the single- spin transitions, 
every term of the density matrix can be determined [|50|]. In 
our case, seven repetitions of the experiment, each with a dif- 
ferent readout pulse, are sufficient to reconstruct the density 
matrix J5l|]. Using this method to determine p out , the projec- 
tions and attenuated correlations, averaged over three different 



input states, and under the action of six different pulses, are 
measured and listed in Table |l[ The input states considered 
are of the form 

Pin=l}+tf+lf, (25) 

for j = x,y,z. 

Because the pulses are short with respect to the decoher- 
ence times, which are all greater than 200ms, the difference 
between the projection and the correlation indicates incoher- 
ent evolution (e.g., RF inhomogeneity) has caused signifi- 
cant errors. This is demonstrated by the fact that the ratio 
of the correlation to the projection, or the attenuation, is on 
the order of the projection (see Table [n|). This further sup- 
ports the need to address errors caused by inhomogeneous ef- 
fects. Robustness to inhomogeneity can be added as a crite- 
rion in pulse determination; the additional resource require- 
ments (time, power, etc.), however, are not known. Decreases 
in the coherent errors are also desirable, but such improve- 
ments must first be sought through improved experimental im- 
plementation. While the goal of fault tolerance clearly hasn't 
been met, these results indicate that experimental implemen- 
tation is nearing ideal simulation results. 



V. CONCLUSIONS 

The ability to implement faithfully a desired unitary trans- 
formation is at the heart of any future implementation of a 
quantum computer. Any control technique should minimize 
the effects of decoherent errors while retaining the required 
addressability. We have demonstrated a method to find con- 
trol sequences that use detailed knowledge of the system's 
Hamiltonian and high-power pulses to create a desired gate by 
strongly modulating the system's dynamics. These gates are 
short in duration yet selective, implementing the correct uni- 
tary transformation. The effect of these control sequences has 
been simulated under various conditions and experimentally 
demonstrated by NMR. Finally, avenues for future improve- 
ments are within the reach of current technologies. 
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TABLE I. Summary of the relevant characteristics for an example 
set of transformations. The three columns list the pulse duration (in 
fis), maximum power (in kHz), and the gate fidelity of the simulated 
pulse. While the maximum power is relatively large, all powers are 
experimentally feasible. The pulses designed for the Crotonic acid 
sample require longer times and yield lower fidelities due to the de- 
creased chemical separation and increase of coupling strengths. 

Pulse Pulse Time(/is) Max. Power(kH z) Fidelity 



Alanine Pulses 



7T] 1 

2 J x 


202 


7.9 


0.9995 


fl 2 

2 J x 


221 


9.3 


0.9995 


f] 3 


212 


9.0 


0.9995 


7T]12 

2\x 


194 


8.5 


QQQS 
yj.yyy j 


tt]23 
2\ x 


179 


9.2 


n none 


2 J x 


163 


10.3 


0.9995 


1 12 


252 


8.0 


0.9996 


^l 23 


129 


10.3 


0.9997 




Crotonic Acid Pulses 




7T] 1 

2\ x 


389 


8.5 


0.9930 


7T]2 

2\x 


610 


4.8 


0.9957 


tt] 3 
2 J x 


392 


6.3 


0.9950 


tt] 4 
2Jx 


559 


7.6 


0.9923 




326 


9.0 


0.9963 




315 


11.3 


0.9932 




345 


8.4 


0.9962 
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TABLE II. Summary of experimental data for the Alanine sample. For each of six tt/2 pulses, the experimentally determined density 
matrix is compared with the expected result. The first column (No Pulse) confirms that the experimental and expected inputs have almost unit 
overlap. Each of the other headings denote which spins are rotated by tt/2. In each case, the projection and attenuated correlation for each 
pulse is averaged over the three inputs ij + if + if, for j = x,y,z. Because all pulses are short in comparison to the natural decoherence 
times, the attenuation gives an indication of the relative significance of the coherent and incoherent errors. Statistical uncertainties arising from 
errors in the tomographic density matrix reconstruction are order 1%. 





No Pulse 


Carbon 1 


Carbon 2 


Carbon 3 


Carbons 1&2 


Carbons 2&3 


All Carbons 


Projection 


0.996 


0.992 


0.995 


0.996 


0.993 


0.996 


0.992 


Correlation 


NA 


0.977 


0.985 


0.986 


0.975 


0.984 


0.986 


Attenuation 


NA 


0.985 


0.990 


0.990 


0.982 


0.988 


0.994 
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C1 


C2 


C3 




7167.8 


54.1 


-1.3 


C1 




-2286.5 


35.0 


C2 






-4881.4 


C3 



Alanine 



Crotonic Acid 



C1 


3767.5 








C2 


72.4 









C3 


-1.3 


69.7 


1915.8 




C4 


7.0 


-1.6 


41.5 


-7850.4 | 




C1 


C2 


C3 



FIG. 1 . Molecular structure and Hamiltonian parameters for both Alanine and Crotonic acid. The chemical shift of each of the carbon nuclei 
is given by the corresponding diagonal elements while the coupling strengths are given by the each of the different off-diagonal elements. 



a) Three-spin pulse 

nT 10 F 1 



b) Four-spin n] pulse 
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2 


3 


Pulse Parameters: 
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4 
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8 
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Time (jis) 


77.3 


78.2 


56.9 


Time (jis) 


44.3 


33.3 


30.1 


16.9 


34.7 


40.9 


37.9 


42.7 


45.4 


Power (kHz) 


8.96 


4.62 


4.06 


Power (kHz) 


9.06 


0.31 


1.36 


3.45 


3.99 


0.94 


1.68 


4.21 


6.42 


Initial Phase (°) 


145 


68 


144 


Initial Phase (°) 


6 


9 


85 


115 


164 


126 


166 


-92 


-34 


Frequency (kHz) 


1.81 


4.59 


9.03 


Frequency (kHz) 


5.66 


0.07 


-5.29 


6.98 


9.30 


-0.04 


0.96 


0.99 


-0.19 



FIG. 2. The ideal RF waveform for two example pulses. The solid (dashed) line is the amplitude (phase) of the waveform. Changes in 
the transmitter frequency (within a single period) were implemented by a discrete linear phase ramp. The sharp discontinuities occur at the 
transitions between periods. Substantial filtering of these high frequency components (smoothing of the shape) has little effect on the gate 
fidelity. In order to experimentally implement the pulse, it is converted into a discrete series of amplitudes and phases (order IK long) by 
sampling the ideal waveform at a constant rate. Details of the pulse parameters (as per Eq. [24]) are listed below each waveform. Due to 
experimental implementation issues, a 6/xs period with zero RF power (i.e., H ex t = 0) is needed before and after the pulse (see waveform) 
and must be included to produce the desired propagator. 
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a) Three-spin7i/2] J pulse 
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FIG. 3. Contour plots of the gate fidelity as each pair of parameters is varied over typical experimental errors, demonstrating that both the 
Alanine and Crotonic acid pulses are most sensitive to errors in the RF amplitude. 



a) Three-spin %/2] 3 pulse 
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b) Four-spin 7i]^ 2 pulse 
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FIG. 4. Gate fidelity of two example pulses as the resonance frequency of a test spin is varied over a range of chemical shifts. The solid 
(dashed) line is calculated with identity (desired transformation) as the theoretical transformation. The vertical dotted lines denote the actual 
chemical shifts for each spin. As can be seen, the gate only works when the test spin is at the appropriate resonance frequency. 
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FIG. 5. Spectra resulting from different sequences of pulses applied to the thermal equilibrium density matrix pthermai = l\ + it + l\. 
All sequences are read from left to right. The reference spectra (resulting from a ^ pulse applied to all 3 carbons) is also shown (grey) for 
reference. Although the chemical shift is order ^ , no significant phase evolution is seen. Selective coupling sequences are also demonstrated. 
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